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I. INTRODUCTION 



QCD evolution equations have been recently established which include pomeron loops 
(or fluctuations) The main feature of these equations, as compared to the Balitsky- 
JIMWLK {3], Q] or Kovchegov equations , are the violation of the geometrical scaling of 
the T-matrix and a modified energy dependence for the saturation momentum {(| . Recently 
the effect of fluctuations also on single diffractive dissociation [7| has been studied and 
compared with previously obtained mean-field-like results |8j at high energies. 

In a recent paper [1] we have studied the QCD pomeron loop equations in zero transverse 
dimensions in which case one can do the calculation of any pomeron loop graph analytically 
(see also The pomeron loop equations are a hierarchy of 



equations which in zero transverse dimensions take the form 

^ = kanW + k{k - l)an (k -V - kf3n^ , (1) 
ay 

where is a normal ordered number operator defined as (n(n — 1) • • • (n — k + 1)) which 
represents the expectation value of /c-pomerons during the evolution. The single terms in 
eq. ([I]) have the following physical meaning: kan^ is the BFKL growth term (a here 
corresponds to ap — 1 = iciaNc In 2 in the real BFKL equation), k{k — l)arS k ~ lS} describes 
fluctuations (pomeron splittings) and kPn^ k+1 ^ recombinations (pomeron mergings). Eq. (CQ) 
is the zero-transverse-dimensional analog of the real QCD equations (see j^). One can easely 
show that ([1]) leads to frame- independent scattering amplitudes when (3 = aa 2 s . 

The scattering amplitudes have been calculated in Ref. l| by treating the recombination 
terms as small perturbations. In a recent paper together with our collaborators 14] we have 
shown that at rapidity Y ^> l/a^a, the perturbative treatment of the recombination terms 
becomes inaccurate and needs to be replaced by more complete calculations (This limit is 
also indicated in our corrections to the LO result given in Eq. (15) in Ref. jl]). Thus, the 
results for diffractive cross sections shown in this paper are only reliable up to Y ~ 1/af. 

In this paper we use the techniques which we have developed in [l] together with the 
AGK cutting rules [17| to calculate single, double and central diffractive cross sections in 
zero transverse dimensions (see Fig. [T]). We do consider only the case where the diffractive 
masses, rapidity gaps and the total rapidity are large. All the relevant pomeron graphs are 
included in the calculations. In four dimensional QCD, the analogous calculations are not 
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yet available. So far only the single diffractive scattering in real QCD including pomeron 
loops has been calculated, however, only in the kinematical region where the diffractive mass 
is not too large [?]. Thus, our calculation in zero transverse dimensions may give indications 
about the diffractive mass dependence of the diffractive cross section in four dimensional 
QCD at very large diffractive masses. We expect the energy dependence of the QCD high 
energy diffractive scattering would bear some resemblances of our result. 






e 



FIG. 1: Typical diffractive dissociation graphs: Single diffractive scattering (A), double diffractive 
scattering (B), central diffractive scattering (C). 

A longstanding problem has been accompanying diffractive dissociation: The lowest or- 
der pomeron contribution to diffractive dissociation disagrees with the Fermilab data (see 



e.g., [181 ]) and it violates unitarity at very high energies. There have been several attempts 
to solve this problem by considering the renormalization of the pomeron flux a larger 
impact parameter or a smaller survival probability 20] at high energies. In this work we 
show that single, double and central diffractive cross sections would naturally fulfill unitarity 
limits when multiple pomeron exchanges are taken into account in addition to the lowest 
order pomeron exchange. 

Furthermore we do find that the differential diffractive cross sections (differential with 
respect to the masses) are suppressed by powers of a s . We will show that this suppres- 
sion is characteristic for differential diffractive cross sections and emerges from pomeron 
splittings/mergings in the t-channel which are needed in order to calculate the differential 
diffractive cross sections. 

The paper is organized as follows: We do start with an illustration of the AGK cutting 
rules, then go over to the single diffractive cross section, which is then followed by double 
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diffractive and central diffractive cross sections. 



II. THE AGK CUTTING RULES 

The AGK cutting rules jl^] (see also 0, represent the generalization of the optical 
theorem for the case of multiple pomeron exchange. Although AGK cutting rules have 
not been proven for QCD, they have been discussed and used in several publications (see 



Refs. 
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24] ) . We will show below how the AGK cutting rules lead to expectable results 
for particle-particle and particle-nucleus results also within this QCD model, and then use 
them to also calculate diffractive scattering processes. 

To illustrate the AGK cutting rules, let us consider for instance the total cross section 
for particle-particle scattering, 



&tot 



2 J2(-l) n - 1 F n , (2) 



?1=1 



where F n is the amplitude for the exchange of n pomerons. According to the AGK cutting 
rules, a tot is related to the elastic cross section <r and inelastic cross section a in , 



Otot = O"0 + Cm , (3) 



which are given by 



= £ F£ (4) 

n=2 



and 1 



oo oo 
k—1 k=l n=k 



£ ** = £ £*?• w 



In the above equations Fq represents the n-pomeron exchange graph where no pomeron 
is cut (elastic cut) as shown in Fig. [2f3 and FJ? is the n-pomeron exchange graph with k 
pomerons being simultaneously cut (inelastic cut) as shown in Fig. [2JA. (cut pomerons are 
marked by crosses). In the next sections we will introduce also the diffractive cut as shown 



1 In eq. (j4]) n starts from 2 because there have to be at least two pomerons, one pomeron on either side of 
the elastic cut. One side of the cut can be viewed as initial state of the scattering while the other side 
can be understood as the final state. 
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FIG. 2: n-pomeron exchange graphs: In graph A A: of the n-pomerons are cut (inelastic cut) and 
in graph B none of the pomerons is cut (elastic cut). 

for example in Fig. [3JA. where only a part of a pomeron is cut. The diffr active cut generates 
diffractive (dissociation) cross sections while the inelastic cut (cut of an entire pomeron) 
leads to a uniform production of final state particles in rapidity. 

For convenience, let us summarize the AGK cutting rules, which allow us to calculate Fq 
and F£ and the diffractive cross sections in the next sections: 

• No matter how many pomerons are cut, there is always one and only one cut which 
indicates the separation between the initial and final states of the scattering; 

• Each cut-pomeron gives an extra factor of (—2) , which can be understood as a result 
of the discontinuity of the pomeron amplitude; 

• Each un-cut pomeron obtains an extra factor of 2 since it can be placed on either side 
of the cut. 

Thus, applying the above AGK cutting rules on Fig. [21 we obtain 




(6) 



and 




(7) 
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where C\ = n\/[k\ (n — k)\] is the number of selections of k pomerons out of n pomerons. It 
is easy to verify Eq. ([3]) by substituting the above expressions for Fq and FJ} in Eq. (jlj) and 
Eq. respectively. 

To gain confidence in the AGK cutting rules, let us show that the cross sections obtained 
using the AGK cutting rules and the pomeron calculus do agree with known unitarity limits 
at high energies. For particle-particle scattering the amplitude for n-pomeron exchanges at 
leading order fl is F» = ,, (the factor (-!)» has aheady been taken into account 

in Eqs. ppp fl)). Inserting this amplitude in Eqs. (l2~]P£f5]) . one finds a tot ~ 2, a ine i ~ 1 and 
o"o — 1 in large Y limit, which is in agreement with the black disk limit. For particle-nucleus 
scattering, now with F n = [a 2 s e aY Vj" , where L is the number of hadrons in the nucleus 
and a 2 L ^> 1, it is straightforward to obtain the following results for the total, elastic 
(diffractive) and inelastic cross sections, 

2a 2 s e aY L 

atot ~ l + a 2 s e« Y V (8) 
2{a 2 s e aY L) 2 

a ° ~ (1 + a 2 e aY L) (1 + 2a 2 e aY L) ' (9) 
2a 2 e aY L 

° m l + 2a 2 s e« Y V l j 

which do as well respect unitarity limits at very high rapidity. 



III. SINGLE DIFFRACTIVE SCATTERING 



The single diffractive process, p + p — > p + X, as sketched in Fig. HJ is a process where one 
of the p-particles breaks up into a " diffractive state" X which is separated by a rapidity gap 
from the p-particle in the final state which remains intact. In the language of pomerons, the 
single diffractive production is generated by a diffractive cut of a multiple pomeron exchange 
graph as shown in Fig. [3A (lowest order graph) and Fig. [3)3 (higher order graph). The 
appropriate variables to describe the single diffractive scattering process are the rapidities 

Y = In ^ and Y = In -^f, where s is the square of the center of mass energy, Mx is the 
diffractive mass, and m p is the rest mass of the particle p. In this picture, the rapidity gap 

Y — Yq is easily understood since there is no cut through any of the pomerons in this rapidity 
window. 

We calculate the single diffractive cross section for the case where Y, Y and Y — Y 
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A 



B 



FIG. 3: Diffractive cut of the multiple pomeron exchange graphs which generate single diffractive 
dissociation: Diagram A is the lowest order graph and diagram B is a higher order graph. 

are large. In this kinematical window the dominant graphs are those shown in Fig. [3 
The leading order contribution of the dominant graphs is calculated in this work. (The 
subdominant graphs which we neglect are those where pomeron loops are formed in the 
diffractive region < y < Yq or in the rapidity gap region Yq < y < Y. The loops in these 
windows are smaller as compared to the dominant loops stretching over the whole rapidity 
Y and give therefore small corrections, proportional to exp (—aY ) or exp [—a (Y — Yq)], in 
comparison to the main result coming from the large loops. In addition there are also other 
subleading contributions neglected here which are a s suppressed as compared to the main 
result which are discussed in detail in Ref. 

According to the topology of Fig. [3] and the AGK cutting rules summarized above, we 
get for the differential single diffractive cross section 2 



where the factor 2 n 1 indicates the number of different diagrams when the diffractive cut is 
made and F n is the uncut single diffractive amplitude for the case of n-pomeron exchange. 

2 We have daso/dYo instead of asD in Eq. (fTTj) since we have considered the production of a particular 
diffractive mass Mx (or Y = \n{M x /mi)). 




(11) 
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We extract F n from the T-matrix for the uncut graphs. The latter is calculated within 
the zero transverse dimensional model [1| which we have recently developed. In the limit of 
exp (aY ) ^> 1 and exp [a (Y — Y Q )] ^> 1, the dominant term is the one where k pomeron 
splittings occur at Y — 0, 

4 fc+1) (F) ~ (A; + l)!e( fc+1 ) ay , 
the last pomeron splitting occurs at Y , 

n[ k+2 \Y,Y ) ~ (* + 2)(* + l) [ Y adyn^ +1 \y)e^ a( - Y ^8(y-Y ) (12) 

Jo 

~ a (k + 1) (Jfc + 2 )\e {k+2)aY ~ aYo , (13) 

and all the k + 1 pomeron mergings occur at Y" 

4 1} (y,y ) - (-l) fc+1 a(fc + l)(A; + 2)!x /c e 2ay - ayo , 

where x = a 2 e aV and T fc (1) (Y,Yo) = a 2 n^(Y,Y ). Now, from T fc (1) (F,F ) = 
(-l) fc ct(£; + 1) (£; + 2)!a Jc+2 e- ay °, we extract F n = a (n - 1) n\x n e~ aY by using F n = 
(— l) n T^} 2 , which inserted into the diffractive cross section given in Eq. (TTTT) yields 

da SD (Y, Y ) ^ ae _ aYo ^ ( _ i)n T _ x {n _ i} (M) 



1 _„,v_ •) <^ 



(15) 



where we have used the definition z = 2a 2 e aY . The first term of the sum in Eq. (|14p 
reproduces the energy dependence of the first order Regge prediction for single diffractive 
processes, and other terms in the summation are contributions of higher order dominant 
graphs. In the high energy limit, when the diffractive mass (or Yq) and the rapidity gap 
(Y — Y ) are kept large, we find 



da SD (Y,Y ) ^ l__ aYo 



dY 2 
or, in terms of the diffractive mass 



—oce 



i 1 + 2 lE - 2 In (2a 2 s e aY ) 



2ale 



2 a aY 



(16) 



dMl 2 \ Ml 

where 7^ = 0.577 is the Euler constant. 



M2 du SD ( S ,Ml) 1 /mgV (17) 



s 



Here, we would like to comment that Eq. (TT6T) is approximately true when olYq and 
ct(Y — Yq) are relatively large, where we find the result ^ae~ aY °. The expectation that 
all inelastic diffraction vanishes for Y — > oo [7] can not be shown in our model since our 
model does not apply for Y ^> 1/af. The origin of our result is obvious: There are (k + 1) 
pomerons in the region < y < Yq. At Yq another pomeron is produced through pomeron 
splitting enforced by the delta-function in Eq. (fl2|) . see also Fig. [3j adding a contribution 
ae a ( Y - Y o) to the (k + 1) pomerons of length Y. This explains the appearance of ae~ aY °. 
The remaining e ay -dependence together with the other pomerons of the same length Y do 
give a sum over all pomeron exchanges which naturally unitarizes the single diffractive cross 
section at large Y. The factor of 1/2 in (Tl6l) is a combinatorial factor emerging after the 
diffractive cut is made. 

Another phenomenon not observed so far elsewhere is that the single diffractive cross 
section in (fl6l) is suppressed by a factor of a. This factor appears through the pomeron 
splitting at Yq. A suppression by powers of a is characteristic for differential diffractive cross 
sections because pomeron splittings in the t-channel (see Figs. I31H1I51) are needed in order to 
calculate them. In the next sections we show how also double and central diffractive cross 
sections are suppressed by powers of a. 

It is tempting to integrate Eq. (fl6l) over Y Q from to Y to obtain a formula for <jsd(X)- 
However, this wouldn't be right since Eq. ffl6l) is not valid for small values of diffractive 
masses and rapidity gaps. When diffractive masses or rapidity gaps are small, many other 
graphs have to be considered besides the dominant graphs shown in Fig. [3j In this case the 
resummation of all order graphs is no longer under control within our formalism. 

In Ref. [3] the (differential) single diffractive cross section in four dimensions including 
pomeron loops has been calculated. However, the calculation is only applicable in the region 
where the diffractive mass Y is small, Y <C 7i/(a s N c ) ln(7r 2 /a^). On the other hand, our 
result applies when the diffractive mass is large, exp[aF ] 3> 1- It would be interesting to 
see whether our result remains valid also in four dimensional QCD. 

IV. DOUBLE DIFFRACTIVE SCATTERING 

The double diffractive scattering process, p+p — > X+X' (see Fig.[Tj3), is generated by the 
three topologically different types of diffractively cut multiple pomeron exchanges as shown 
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FIG. 4: Three topologically different types of the diffractively cut multiple pomeron exchange 
graphs which generate double diffractive dissociation. Diagram I (type I) is the lowest order graph; 
Diagram II (type II) is a higher order graph which requires at least 3 pomerons exchanged; Diagram 
III (type III) is another higher order graph which needs at least 4 pomerons exchanged. 

r-n M f ^ 

in the lowest order in Fig.Hl We use below the variables Yq = In —4- and Y\ = Y —Y[ = In —4-, 
where Mx and M' x are the diffractive masses, to express the double diffractive cross section. 
The rapidity gap between the diffractive states X and X' is Y — Y Q — Y\. 

Let us first compute the dominant contribution of the un-cut amplitude for the case of 
multiple pomeron exchange with one pomeron splitting at Yq and one pomeron merging at 
Y\. Further we consider the case where Yq and Y\ are kept large. The result after the last 
pomeron splitting at Yq which is given by Eq. (IT2|) . when followed by a pomeron merging at 
Y{, gives 

nt +1) {Y{, Y , Y) ~ -a] (k + 1) f adynt +2 \y)e( k+1 W Y ^6 (y - Y{) , (18) 

Jo 

= (-a 2 s ) (k + 2)! (k + if a ^+^Y+aY{-aY ^ (1Q) 

and after k more pomeron mergings at rapidity Y (see fl|), one ends up with Tj^\Y, Yi, Y ) ~ 
(-l) fc+1 a 2 (k + l) 2 (k + 2)\x k+2 e- aYl - aYo , where Y 1 = Y — Y{. 

The resulting F^ D = a 2 in — l) 2 n\x n e~ aY °~ aYl has the following interpretation: There 
are (n — l) 2 different graphs in the n-pomeron exchange case, each of them contributing 
by the same amplitude a 2 n\ x n e _ay °~ QYl . Among these (n — l) 2 pomeron exchange graphs, 
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there are three topologically different types when the diffractive cut is made (e.g., see Fig.H]). 
In the n-pomeron exchange case, it is easy to find that there are = (n — 1) diagrams 
of the first type, Njp = 2 (n — 2) diagrams of the second type and Njfj = (n — 2) (n — 3) 
diagrams of the third type. (Evidently, + + Njjj = (n — l) 2 .) When the cut 
is made, the number of diagrams of type I becomes N\ n ^ ut = (n — 1) 2 n_1 , the number of 
diagrams of type II turns out to be Njj cut = (n — 2) 2 n ~ 2 and the number of diagrams of 
type III is Njjj cut = (n — 2) (n — 3) 2 n_3 . With the total number of diagrams after the cut, 
Ncut = [4 (n — 1) + (n — 1) (n — 2)] 2 n_3 , we get for the double diffractive cross section 

da DD (Y, Y , yQ _ 2 Yo _ aYl ^ n n {n ) 
dY dY 1 ~ ^ ' 

In high energy limit, while Y and Y\ large, one obtains 

da DD (Y, Y , yQ _ I 2 Yo _ aYl 
dY dY 1 ~ 4 



or, equivalent ly, 



M 2 M l2 daDD ^ M ^ M '^ ~ la? f "iY f ™i 



dMldM% 4 \Ml) \M'^ 

The double diffractive cross section equals the product of two single diffractive cross 
sections at very large Y, which is in line with the naive expectation dapDj 1 dY§dY\ ~ 
[dasD/dY dosol dY^\j o t ot since the double diffractive process can be viewed as two sep- 
arate and independent single diffractive processes occurring at both, the projectile and the 
target. Therefore one obtains a stronger a and diffractive mass suppression for the differen- 
tial double diffractive cross section as compared with the differential single diffractive cross 
section. Also the double diffractive cross section unitarizes due to the multiple pomeron 
exchanges, instead of showing a rapid growth, when Y becomes very large large. 



V. CENTRAL DIFFRACTIVE SCATTERING 



The central diffractive scattering, p + p — > p + X + p (see Fig. (DC), as shown in Fig. [5] 
(higher order graph), is obtained by following the same procedure as in the previous sections. 
The un-cut amplitude of k-th order (i.e., k pomeron exchanges between and Y in addition 
to the lowest order diagram) reads 

Tf } (F, Yl, Y ) ~ (-l) k+2 2a 2 a] (k + 1) (k + 2)\x k+2 e - aYl+aY ° . 
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Y 



FIG. 5: Diffractive cut of general multi-pomeron exchange graphs which generate central diffractive 
dissociation. 

When a diffractive cut is made, the diffractive cross section becomes 



and for high energy scattering but Y\ — Y Q = 5Y = In y-^rj large (Mx is the central 
diffractive mass) it reduces to 



The suppression ot 2 a 2 s can be easily understood by looking at the lowest order graph (Fig. [5] 
without the two pomerons on both sides of length Y): The merging at Yq gives an (3 and the 
subsequent splitting at Yj. a a. The rapidity dependence is exp[aY] exp[a(Y— (Yi— Yq))]. The 
remaining factor of comes from the coupling of the pomerons to the scattering particles, 
thus, giving altogether a 2 a 2 exp[— a{Y\ ~ ^o)] x<2 , where x = a 2 exp[aY], which explains the 
result in fl20l) . The higher order graphs obtained by adding further pomeron exchanges as 
shown in Fig. [5] lead to a Y-dependence in accordance with unitarity limits. The differential 
central diffractive cross section is strongly a s -suppressed (a ~ a s ) as compared with the 
differential single (double) diffractive cross section. 

In summary, we systematically calculate single diffractive, double diffractive and central 
diffractive cross sections within the zero transverse dimensional Regge model [lj. The resum- 






(20) 
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mation of all higher order dominant diagrams shows that diffractive cross sections gradually 
approach unitarity limits instead of increasing rapidly at very large energies. Moreover, we 
do find powers of a suppression in the inclusive diffractive cross sections. 

Acknowledgements: We acknowledge numerous stimulating discussions with Alfred 
Mueller. A. Sh. acknowledges financial support by the Deutsche Forschungsgemeinschaft 
under contract Sh 92/2-1 and wishes to thank the theory group at DESY for hospitality 
during his visit when this work was being completed. 



[1] A. I. Shoshi and B. W. Xiao, Phys. Rev. D 73 094014 (2006). 

[2] A. H. Mueller, A. I. Shoshi and S. M. H. Wong, Nucl. Phys. B 715 440 (2005); E. Iancu and 
D. N. Triantafyllopoulos, Phys. Lett. B 610 253 (2005); A. Kovner and M. Lublinsky, Phys. 
Rev. D 71 085004 (2005); Y. Hatta, E. Iancu, L. McLerran, A. Stasto and D. N. Triantafyl- 
lopoulos, Nucl. Phys. A 764 423 (2006). 

[3] I. Balitsky, Nucl. Phys. B 463 99 (1996); Phys. Rev. Lett. 81 2024 (1998); Phys. Lett. B 518 
235 (2001). 

[4] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. Phys. B 504 415 (1997); 

Phys. Rev. D 59 014014 (1999). E. Iancu, A. Leonidov and L. D. McLerran, Phys. Lett. B 

510 133 (2001); Nucl. Phys. A 692 583 (2001); H. Weigert, Nucl. Phys. A 703 823 (2002). 
[5] Y. V. Kovchegov, Phys. Rev. D 60 034008 (1999); Phys. Rev. D 61 074018 (2000). 
[6] A. H. Mueller and A. I. Shoshi, Nucl. Phys. B 692 175 (2004); E. Iancu, A. H. Mueller and 

S. Munier, Phys. Lett. B 606 342 (2005). 
[7] Y. Hatta, E. Iancu, C. Marquet, G. Soyez and D. N. Triantafyllopoulos, "Diffusive scaling and 

the high-energy limit of deep inelastic scattering in QCD at large N(c)," Nucl. Phys. A 773, 

95 (2006). 

[8] S. Munier and A. Shoshi, Phys. Rev. D 69 074022 (2004). 

[9] A. H. Mueller and G. P. Salam, Nucl. Phys. B 475 293 (1996); 
[10] P. Rembiesa and A. M. Stasto, Nucl. Phys. B 725 251 (2005). 
[11] M. Kozlov and E. Levin, Nucl. Phys. A 779, 142 (2006). 

[12] M. Kozlov, E. Levin, V. Khachatryan and J. Miller, "The BFKL pomeron calculus in zero 
transverse dimensions: Diffractive processes and survival probability for central diffractive 



13 



production," |arXiv:hep-ph/0610084 



[13] J. P. Blaizot, E. Iancu and D. N. Triantafyllopoulos, "A zero-dimensional model for high- 
energy scattering in QCD," |arXiv:hep-ph /0606253, 

[14] S. Bondarenko, L. Motyka, A. H. Mueller, A. I. Shoshi and B. W. Xiao, "On the equiva- 
lence of Reggeon field theory in zero transverse dimensions and reaction-diffusion processes," 



arXiv:hep-ph/0609213, 



[15] A. Kovner and M. Lublinsky, Nucl. Phys. A 767 171 (2006). 

[16] E. Iancu and D. N. Triantafyllopoulos, Nucl. Phys. A 756, 419 (2005). 

[17] V. A. Abramovsky, V. N. Gribov and O. V. Kancheli, Yad. Fiz. 18 (1973) 595 [Sov. J. Nucl. 

Phys. 18 308 (1974)]. 
[18] K. Goulianos, Phys. Lett. B 358(1995) 379 . 

[19] A. H. Mueller, Eur. Phys. J. A 1 19 (1998);A. Hebecker, Phys. Rept. 331 1 (2000). 
[20] E. Gotsman, E. Levin and U. Maor, Phys. Lett. B 438 229 (1998). 

[21] J. Koplik and A. H. Mueller, Phys. Rev. D 12 3638 (1975); E. Levin, "An introduction to 

p omer ons , " arXiv : hep-ph /9808486. 
[22] J. Bartels, M. Salvadore and G. P. Vacca, Eur. Phys. J. C 42, 53 (2005); J. Bartels and 

M. G. Ryskin, Z. Phys. C 76, 241 (1997). 
[23] Y. V. Kovchegov and E. Levin, Nucl. Phys. B 577, 221 (2000); M. A. Braun, Phys. Lett. B 

483, 105 (2000). 
[24] D. Treleani, Int. J. Mod. Phys. A 11, 613 (1996). 



14 



